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1. Introduction 

In the 35 years since Hawking's seminal paper on the quantum emission of radiation 
from a black holepQ, a large body of work has been devoted to solving the so called 
information loss problem. For a black hole of fixed mass M, this emission process 
yields a black body spectrum with characteristic temperature Th = hc 3 /8nkBGM, 
irrespective of the initial state of the matter from which the black hole is formed. The 
inability to reconstruct the initial, possibly pure state of the black hole from the total 
emitted radiation signals the apparent breakdown of unitary evolution and the ^-matrix 
description of the Hawking process. Although Hawking's calculations have since been 
verified in a number of ways[2j El HJ El E] , this breakdown at the foundation of quantum 
mechanics suggests that our understanding of black hole dynamics is not yet complete. 

The information loss problem rests on two key assumptions made in the standard 
derivation of Hawking radiation: (1) the perfectly thermal (i.e. mixed) character of 
the outgoing radiation; (2) the validity of this emission process over the lifetime of 
the black hole. The notion of thermal spectrum considered here is not that of a 
blackbody frequency spectrum, but rather the quantum thermal probability distribution 
defined by the temperature Th of a single mode (single frequency) of Hawking radiation. 
The traditional picture of the Hawking process leaves no room for deviations from 
this thermal distribution and thus breaks the requirement of pure-state— >-pure-state 
evolution enforced by unitarity. By itself, this process need not lead to information loss 
since the information content of the black hole may be stored in entanglement between 
particle pairs created on opposite sides of the horizon[7J [8]. Results from many-body 
theory suggesting that entanglement across a boundary scales with the area of the 
boundary lends credence to this view|9J. However, with the second assumption, the 
black hole causes a loss of entanglement and thus information. 

With the expectation that information must be conserved, many suggestions for 
resolving the information loss problem have been put forward. Current proposed 
solutions include long-lived and stable Planck-scale remnants [Ml EI], baby universes|12[ 
[T3] . and the possibility of information escaping as non-thermal Hawking radiation 
[TH [To] . In all of these proposals, corrections to the Hawking process manifest 
themselves when the black hole has evaporated to a size near that of the Planck length, 
l p = ^JhG/c 3 ~ 10~ 35 m, at which quantum gravitational effects, neglected in Hawking's 
original analysis, are expected to play a role. In considering quantum states of the 
gravitational field, there is the possibility of back-reaction and entanglement of the 
radiating matter degrees of freedom with those of gravity [16, 17J. Although it is natural 
to consider a quantized gravitational field for the Hawking process, the current lack of 
a full quantum mechanical description of gravity severely limits progress in directly 
addressing this scenario. In fact, exactly which degrees of freedom, if any, should 
be quantized is still subject to debate [TH [19]. In this paper we investigate a simple, 
zero-dimensional quantum optics model of Hawking radiation that mimics some of the 
essential physics present in the original information loss problem. 
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As a zero- dimensional model we consider the following trilinear Hamiltonian: 

H = hu a a + a + hubb + b + hu c c + c + i%x (ab + d + — a + b&j , (1) 

consisting of three harmonic oscillator modes with the frequency relation, u a = Ub + oj c . 
We will designate the modes as pump (a), signal (b), and idler (c) respectively. This 
Hamiltonian describes several quantum optics processes including frequency conversion, 
Raman and Brillouin scattering, the interaction of two-level atoms with a single 
mode resonant EM field, and is the full quantum generalization of the parametric 
amplifier[20l ED E21 E31 122 E3 E51 EZJ. The connection between black hole radiance 
and parametric amplification was appreciated shortly after Hawking's discovery [28J. 
Both processes amplify vacuum fluctuations resulting in the production of correlated 
photon pairs. Tracing over one of the two subsystems (i.e., signal and idler) yields 
statistics that are identical to that of a thermal distribution [29, 30J. The energy 
source in the parametric amplifier is assumed to be a classical pump such as a laser 
or microwave generator with fixed amplitude driving a system with nonlinearity, 
the first nonlinear susceptibility in a medium without inversion symmetry |31j. Viewed 
as a black hole model, the pump plays the role of black hole mass M while the signal 
and idler modes of the parametric amplifier represent the escaping and trapped Hawking 
photons, respectively, as depicted schematically in Fig. [T] The trilinear Hamiltonian (IT]) 




black hole(a) 
(a) (b) 



Figure 1. a.)Diagram depicting the dynamics of the trilinear Hamiltonian. Initial 
vacuum modes are not shown for clarity. The nonlinear interaction is generated by a 
system possessing a second-order susceptibility x*- 2 '. b.)Equivalent dynamical elements 
involved in the Hawking process. 

generalizes the parametric amplifier by quantizing the pump mode and allowing for 
energy loss to the signal and idler modes. The expectation value of the pump mode 
energy is analogous to the mass M of a quantum mechanical black hole. We will explore 
this model by establishing the conditions under which the signal mode spectrum of the 
trilinear Hamiltonian deviates from the predicted thermal spectrum of the conventional 
parametric amplifier. We shall see that the quantization of the pump mode degree of 
freedom results over time in entanglement with the signal and idler modes and dynamics 
that become markedly different from the parametric approximation. In particular, the 
signal mode develops a strongly non-thermal spectrum that is dependent on the initial 
pump mode state. The corresponding entropy is reduced relative to that of a thermal 
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(maximally mixed) state indicating the presence of information. These model system 
results lend support to the view that late-time Hawking radiation contains an increasing 
amount of information about the initial quantum state of the black hole and is composed 
of particles entangled with quantized gravitational states. 

The outline of this paper is as follows: In Sec. [2] we review the derivation of 
the amplification of vacuum fluctuations under the parametric approximation. Sec. [3] 
considers the semi-classical approximation whereby the back-reaction from the quantized 
radiation onto the classical pump is accounted for and derives the self-consistent 
equations of motion for the signal and idler modes. In Sec. |4]we consider the full quantum 
dynamics Eq. ([TJ under the short time analytical approximation and compare this result 
along with the full quantum numerical solution to the parametric and semi-classical 
approximations from the previous sections. Sec. [5] investigates the role of entropy and 
entanglement in the production of non-thermal states. Finally, Sec. [6] concludes with a 
brief discussion of the results and consequences for black hole evaporation. 



2. Parametric Amplifier and Hawking Emission 

In the following, we derive the well-known thermal spectrum of the signal mode under 
the parametric assumption of a fixed amplitude pump mode. Replacing the pump mode 
in Eq. (jl]) with a fixed amplitude drive A results in the interaction frame Hamiltonian 

Hi = ih x A (b + c + - be) . (2) 

The Heisenberg equations of motion for the signal and idler mode operators are, 

d -f = A X cm d -f = A xm +. (3) 
These can be readily solved yielding the Bogoliubov transformations, 
b( T ) = 6(0) cosh (At) + c(0) + sinh (At) 

( 4 ) 

c (r) = c(0) cosh (At) + b{0) + sinh (At) , 

where we have expressed the dynamics in dimensionless time r = xt- If the system 
starts with both signal and idler modes in the ground state, |^(0)) = |0, 0) bc , then 
Eq. (El) gives for the number operators Nb and N c : 

N b (T) = N c (t) = sinh 2 (At) . (5) 

Additionally, we are interested in the probability distribution of the individual signal 
and idler subsystems in the number state basis. With the system initially in the ground 
state, the unitary evolution corresponding to Eq. ^ can be expressed as, 

|*(r)) = exp \At (b + c + - be)} |0, 0) bc . (6) 



Making use of the disentangling theorem [32], 

exp \At (b + c + - be)} = e ™+c+ e - 9 ( b +6 +c+c+ i) e -r bC; (?) 
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with T = tanh (At) and g = In cosh (At), where the last term in Eq. ^ vanishes for 
the ground state and the middle term reduces to e~ 9 , we have 

T n (b + c 

— I u > u /6c = 

n=0 



|o,o) 6c = E 



n=0 



TV. 



|0,0) 6c = ^tanh™ (At) \n,n) 



be ' 



(8) 



with the evolving state vector \^(t)) given by 

|\J(r)) = sech (At) ^ tanh" (At) \n, n) bc . 



(9) 



Let us now focus on operators acting on a subsystem spanned by the states of either the 
signal or idler mode individually. Assuming we are interested in the signal mode only, 
tracing over the idler subsystem in Eq. ^ gives the operator expectation value 

(0{r)) b = sech 2 (At) £ tanh 2n (At) (n\0\n) h (10) 



for an arbitrary signal mode operator O. Comparing Eq. (10) with the spectrum of a 
thermal state defined by temperature T, 

(O) = Y,Pn (n b \0\n b ) = (l - e - h ^ k » T ) £ e ~ nh ^ T (n b \0\n b ) , (11) 

n n 

indicates that the signal mode is in a thermal state provided we define the temperature 

as 

T ^ = 2k B In [coth (At)] ' (12) 

where the time dependence is a consequence of the rapid increase in the occupation 
number of the modes, Eq.Q. As in the Hawking process[33], the particle pairs 
generated by the parametric amplifier form an entangled two-mode squeezed state given 
by Eq. (|9| [34] . The bipartite structure of this system allows calculating the entanglement 
between particle pairs via the Von- Neumann entropy Si, also referred to as entanglement 
entropy. For mode i — b,c, using the reduced density matrix pi we have 

Si = -Tr (pi In , (13) 

where we drop the usual k-Q factor. For both the signal and idler modes, this entropy is 
given by 



S 



th 



In 



-hu/k B T(T) 



1 



hw/k B T(T) 



-1 



(14) 



k B T(T) 

which is the thermal entropy of a quantum harmonic oscillator with temperature defined 
by Eq. (12). Thus, as for the Hawking thermal radiation, we see that the temperature 
of the parametric oscillator is determined by the entropy generated from tracing over 
one of the two modes in a particle pair squeezed state. 
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3. Semi-Classical Analysis: Backreaction on a Classical Pump Mode 

We now go beyond the just-considered parametric amplifier model and incorporate 
backreaction effects due to the emission process using a semi-classical approximation 
where the pump mode is treated as a classical variable, while the signal and idler modes 
are quantized. In order to self consistently solve for the system evolution, we first work 
out the dynamics of the pump mode assuming it behaves as a classical variable affected 
by the expectation values of the signal and idler modes, after which we substitute 
the resulting c-number expressions for the pump operators into the Hamiltonian and 
calculate the signal and idler evolution. A similar procedure arises for the semi-classical 
Einstein equations: 

^ = ^(^1*), (15) 

where the mass M black hole spacetime geometry is modeled classically through the 
usual Einstein tensor interacting with the quantum matter/radiation fields via 

the expectation value of the stress-energy-momentum tensor operator T^ v with respect 
to a suitable incoming quantum field state Solving these semiclassical equations 
using an analogous procedure to that outlined above results in a nonlinear Schrodinger 
equation [35] and also a possible loss of spacetime stability [36]. In addition, issues such 
as non-causal dynamics [37] and possible inconsistencies with wavefunction collapse [38J 
occur. In contrast, the semi-classical model presented here can be readily solved by 
exploiting the symmetries present in the Hamiltonian. 

To begin, we consider Eq. ([I]) in the interaction frame, 

Hi = ih X (ab + c + - a + bc) (16) 

and obtain the following mode equations, 

da db + dc + 

= - x bc; — = xac + ; — = xab + (17) 
at dt dt 

that lead to the evolution of the pump mode number operator, 
dN a 
dt 

It is easy to show the corresponding number operators for signal and idler modes are 
given by dN b /dt = dN c /dt = —dN a /dt. To proceed, we will use the following Manley- 
Rowe constants of motion 



X (ab + c + + a + bcj . (18) 



M ab = N a + N b - M ac = N a + N c ; M bc = N b - N c , (19) 
expressing the underlying SU(2) and SU(1,1) symmetries in our model 001 HI]- 



Differentiating Eqs. (17) again results in decoupled equations of motion containing only 
commuting operators. 

d 2 N r 

a - 2 X 2 \3N 2 a - N a (2M ab + 2M ac + 1) + M ab M ac 



dt 2 
d 2 N b 

d 2 N c 



-2 X 2 [3N 2 - N b (AM ab - 2M ac - 1) + M ab (M ab - M ac - 1)] (20) 
-2 X 2 \3N 2 - N c (4M ac - 2M ab - 1) + M ac (M ac - M ab - 1) 
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We now take the expectation value of the pump number operator equation (20) for N a 
and make the approximation (N%) « (N a ) (N a ), the validity of which can be checked 



using full quantum numerical simulations of Eq. (17) [see Sec. [4]. These approximations, 
along with the condition that both signal and idler mode start in the vacuum state, lead 
to the semiclassical evolution for the pump mode: 



N a (r) = (3 + + [N a (0) - f3 + ] dn 



" J-r., 



iVa(O) ~ ft. 



where dn(u, k) is the Jacobi elliptic function and 
1 



ft 



± 



1 + 2N a (0) ±Jl + 12N a (0) + 4N a (0y 



(21) 



(22) 



It is important to note that both Eqs. (21) and (22) are expressed in terms of the initial 
conditions of the pump mode only, a consequence of the relations (19) and the fact that 
the signal and idler are initially in the vacuum (ground) state. Equations of motion for 
both signal and idler can then be obtained by substitution of the c-number expressions 
for the pump mode with amplitude given by Eq. (21), 



N a (r)e 



i<j>(t) 



(23) 



where <j>{t) is a slowly varying function of time, into Eq. (16): 

H l = ihxNait) 1 ' 2 (6+c+e-^W - bce^V) . 

The time evolution of this now bilinear Hamiltonian can be straightforwardly calculated 
if we assume 



(24) 



N a (t) 



d(j>{t) 
dt 



const, 



(25) 



which is the most general condition under which the Hamiltonian at different times 



commutes, 
relation 



, da da + 

a — a— — 

dt dt 



0. Combining Eqs. |l7j ) and (23) results in the phase 

d(j){t) 



—2iN„ 



dt 



(26) 



which, expressing the left-hand side as Eq. (16) through the Heisenberg equations for 
the operators Eq. (17), yields a second condition for the pump phase: 



2N a (t) 



d(j)(t) 
dt 



(27) 



As long as the initial state contains at least one mode in the vacuum state, we have the 
constant of motion (Hi) = 0. Thus, Eqs. (25) and (27) are only consistant if <f>(t) = 0. 
This leads to a bilinear Hamiltonian of the form, 

H T = ih X N a (t) 1/2 (b + c + - be) . (28) 

Using a similar analysis to that employed in Sec. [2] gives, 
db{r) 



dr 
dc{r) 

dr 



b(0) cosh (0(r)) + c + (0) sinh (0(r)) 
c(0) cosh (0(r)) + 6+(0) sinh (0(r)) 



(29) 
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where, 

6(t) = jT ^N a (r')dr'. (30) 

Therefore, we see that each mode behaves similarly as in the parametric approximation 
but with the number of particles dependent on the pump dynamics: 

N b (r) = N c (t) = sinh 2 (0(r)) . (31) 

However, most importantly the modes are again in a thermal state with distribution 

(0(r)> = sech 2 (0V o (t)t) X>nh 2n (0V^)t) (n|0|n) . (32) 

Thus, the spectrum of the signal mode (Hawking radiation) remains thermal when pump 
mode energy loss due to signal and idler particle creation is taken into account within 
the semi-classical approximation. 

4. Pull Quantum Description 

4-1- Short-time Approximation 

In this section we consider the full quantum dynamics of the trilinear Hamiltonian 
Eq. (fl|, (i.e. quantize the pump mode as well). Although the trilinear system does not 
allow for a general analytical solution, we can obtain approximate equations in the short 
time limit r = %i <C 1. Using the condition u b = u c = u a /2, appropriate for modeling 
particle generation by a black hole, we can rewrite Eq. ignoring constant terms as 

H = H + H int = hco a (a + a + K z ) + ih X (aK+ - a + K_) , (33) 

where 

K + = b + c + , K- = be, K z = ^ (b + b + c+c + l) . (34) 
We now construct the Casimir operator, 

K 2 = K 2 Z - 1 - (K + K^ + K_K + ) = 1 (Ml - l) , (35) 



where the second equality comes from definition ( 19 ). This operator obeys the eigenvalue 
equation 

K 2 \m) = k(k-l) I*), (36) 

where k = 1/2 (|Mj, c | + 1). We now look for simultaneous eigenstates of both K 2 and K z : 
\k,n c ), where n c is the number of idler particles. These states can also be decomposed 
using the number state basis 

\k,n c ) = \n c + 2k- l) b \n c ) c = \n b ) \n c ) , (37) 

with nfy representing the number of particles in the signal. Operating on this state with 
K z gives 

K z \k, n c ) = (k + n c ) \k, n c ) = ^ — + ^ c + \n b ) \n c ) , (38) 
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valid only for initial conditions where both signal and idler modes are in pure states 
satisfying n b > n c . Including the pump mode state in the full state vector and expressing 
the idler mode population as a function of the pump amplitude 

\n a )\k,M ac -n a ) , (39) 



one may switch to the interaction frame where the evolution depends only on H int (33) 
\t; k, M ac ) = e -K+- a+ *-) |0; k, M ac ) , (40) 

with the initial state with M ac = N a (0) is given by 

\0;k,M ac ) = \M ac )\k,0). (41) 



The short time approximation r = yt <C 1 to Eq. (40) can be calculated using the 
Baker-Campbell-Hausdorff formula truncated to O (t 2 ) 

e r(aK+-a+K„) _ e raK+ e ~ra+ K. g - ^-[aK+,a+ KJ\ _ ^) 

Acting with the third exponential operator term on our initial state gives 

e -4[aK + ,a+K-} ^ ~ A _ kM ac T 2 ) \M ac ) \k,0) (43) 



indicating clearly that the time over which the approximation is valid t = rj% < 
l/(xVkM ac ) decreases with pump amplitude and coupling strength. Evaluating the 



BCH approximated expression (42) on the initial state (41), we obtain for the full 
evolution of the state: 

1 M ac 

\t; k, M ac ) = J2 fn (k, M ac ) r n \M ac -n)\k,n), (44) 

\/ N M ac (r) n=0 

where iVjvf oc (r) is the time-dependent normalization factor 

N Mac (r) = e ^\ 2M -T (M ac + 1, r- 2 ) , (45) 
with T (a, b) the reduced gamma function and 



fn(k,M a 



MJT{2k + n) 12 



(46) 



n\ (M ac -n)\T(2k)_ 

Our interest in the crossover from classical to quantum dynamics for the pump mode 



suggests that we use coherent states built from linear combinations of Eq. (44) for the 
initial pump state. To this end, we consider a general initial state 



CO 



l*(o)> = E^4=l s >l°>l°>> ( 47 ) 



s=0 



N s (0) 



with the pump mode in an as yet unspecified initial pure state with probability P s = \a s \ 
of being in state s and both signal and idler modes in the vacuum state (k = 1/2). Here 
we implicitly assume the probabilities are normalized and add to unity. The density 



matrix at some later time r resulting from (47) is given by 



oo oo s r 



pUt) = E E E E 4==4=S I s *> < r - '"I & ' I ' ( 48 ) 

s =0r=0i=0j=0 J N S (T) J N r (T) 
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which, after performing a partial trace, leads to the reduced density operators 

pt(r) =ttit ^^4=41^^ |s - i) (r - jl ' (49) 

s=0 r=0 i=0 j=o ^/iV s (r) JN r [T) 
where <5j.j is the Kronecker delta, and 

Pf(T)=±±P s P^T 2i \i)(i\ (50) 

for the pump and signal modes respectively. Fig. [2] compares the expectation values of 
the pump and signal modes for the parametric, semi-classical, and short-time quantum 
approximations, as well as the full numerical solution to Eqs. (17) where the pump mode 
is initially in a coherent state with amplitude (iV a (0)) = 9 corresponding to classical 
pump modes with amplitude A = 3, and the signal and idler modes are initially in their 
vacuum (ground) states. 




Figure 2. a.)Evolution of pump mode initially in coherent state (N a (0)) = 9 for the 
parametric (red) , semi-classical (blue) , short-time quantum approximation (green) , and 
full quantum numerical solution (black), as a function of dimensionless time r = 
b.) Corresponding population of the signal mode. Note that the semi-classical analysis 
ceases to be valid when the pump mode becomes depleted. 



4-2. Mode Spectrum Dynamics Under the Short-time Approximation 

The main limitation of the short-time approximation is the inability to account for 
backreaction effects resulting from the build up of quanta in the signal and idler modes. 
For particles produced in the Hawking process, however, the entangled pairs generated 
early in the evolution do not, to first-order, effect those created later from the black 
hole[45j. Furthermore, the emitted radiation does not build up in the vicinity of the 
black hole, but escapes to spatial infinity. Therefore, in this section we will suppose that 



the expressions derived above for the evolving pump and signal states, Eqs. (49) and 



(50) respectively, in fact provide a more relevant zero- dimensional model of a radiating 



black-hole when extrapolated beyond their original short-time domain of validity. The 
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pump and signal modes contain in general a large number of coefficients for each number 
state and cannot be easily evaluated. However, in the long-time limit where the pump is 
nearly depleted, these equations can be considerably simplified by noting that a general 
element of Eq. ( 50 ) for fixed s is given by 

- 2 r( s + i) 



m 

! NJt 



-T 



2 1 



UlT 



2u 



As such, only those elements where i 
density matrix 



Pb 



u = s-i. (51) 

[S+1,T l ) 

= s remain nonzero, leading to the signal mode 



(52) 



which is in general a mixed state with number-state probability distribution P s identical 
to that of the initial pump state. Therefore, by measuring the late-time signal or idler 
mode distribution we recover all but the phase information associated with the initial 
pure state of the pump mode. Additionally, Eq. (52) shows that the number of quanta 
in the signal mode is equal to the initial number of quanta in the pump mode, indicating 
that the pump (49) ends up in the vacuum (ground) state. 



The most important knowledge gained from Eq. (52) is that the signal mode 
spectrum will no longer be that of a thermal state, in contrast to the parametric and 
semi-classical approximations in Sec. [2] and Sec. |3j respectively. Focusing on coherent 
states, in Fig. [3] we give an example of the evolution of Eqs. (49) and (50) by plotting 
the probability distributions for both modes as a function of time r for a pump initially 
in a coherent state with amplitude (N a (0)) = 9 and initial vacuum state for both signal 
and idler. Additionally, in Fig. ^jp we highlight what the thermal distribution would be 




Figure 3. a) Spectrum of the pump mode at several time-steps when the pump is 
initially in a coherent state with (N a (0)) — 9. b) Corresponding spectrum of the signal 
mode. Red lines indicate what the distribution would be at each time-step if in a 
thermal state corresponding to occupation number (Nb(r)), Eq. (53). 



at each time step by equating (N b (r)) with a Bose-Einstein distribution to extract an 
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effective temperature: 
(N b (r)} 



D hu) b /k h T cf{ (T) 



- 1 



(53) 



In Fig. [3^, see we the evolution of the initial coherent state as it loses quanta to the signal 
and idler modes and evolves towards the ground state represented by a peak in the 
probability distribution at the n = number state value. The corresponding evolution 
of the signal mode in Fig. [3)3 starts with the vacuum state and progresses towards the 
state with probability distribution identical to that of the initial pump coherent state, 



Eq. (52). By comparison with the effective thermal state (53), we can see that the 



initial probability distribution for the signal mode is nearly that of a thermal state 
up until the pump mode has transferred nearly half of its initial energy corresponding 
to (Nb(r)) = 4.5. As the evolution continues, Fig. |3|d shows the increasing deviation 
from the effective thermal description for the signal mode distribution as expected from 



Eq. (52) 



4-3. Non-thermal Spectra and Information 



We now quantify the deviations of the signal mode spectrum Eq. (50) from that of a 
thermal state using the fidelity [42J 



F b (r) = Tr v /p 6 (r)V2 <r ( r ) p6 ( r )i/2 > (54) 
where Pb{^) is the density matrix of the signal mode and o~(t) is a thermal density 



matrix with effective temperature determined by Eq. (53) using the occupation number 
(iV 6 (r)). The fidelity provides a measure of the distance between these two states in 
distribution space with range < F < 1, where unity indicates that the two density 
matrices are identical. In Fig. |4^l we plot the fidelity of the signal mode assuming a pump 
mode initially in a coherent state starting with several occupation numbers ranging 
from one to nine. By comparison with Fig. [2j we can see that for the initial coherent 
state (N a (0)) = 9 the fidelity remains essentially unity, indicating that the signal mode 



density matrix Eq. (50) is nearly thermal, until the pump has transferred close to half 



of its initial energy into the signal and idler modes beyond which point there is a strong 
deviation from the thermal state as the signal mode asymptotically approaches the state 



given by Eq. (52). This is in agreement with the qualitative description presented in 



Fig. [3]d and remains true for all the initial states considered in Fig. [4^,. 

This deviation from a thermal distribution also indicates that the signal mode 
spectrum contains information defined as [14] 

I b (r) = St h (r) - S b (r), (55) 

where S^ h (r) is the Von Neumann entropy of the signal mode in a thermal state with 



equal average occupation number, Eq. (14), and »S&(t) is the actual entropy of the mode 



calculated using Eq. (13) and the reduced density matrix of the signal mode p&(r). 



In Fig. 4b we plot the information contained in the signal mode for the initial pump 
coherent states considered in Fig. 4a. Given that the fidelity is nearly unity up until the 
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Figure 4. a) Fidelities of the signal mode given a pump mode initially in a coherent 
state with occupation number (N a ) — 1 (blue), 3 (red), 6 (green), and 9 (black) under 
the short-time approximation, b) Information present in the signal spectrum calculated 



from Eq. (55 1. Equivalently colored vertical lines denote the time at which the effective 

jeff 
1 bc ■ 



subspace dimensions satisfy, c?° ff = df f 



pump mode transfers half of its original quanta to the signal and idler modes, it follows 
that the information content of the signal or idler mode is close to zero over this same 
time interval before increasing as the signal spectrum becomes identical to that of the 
initial coherent states. 

In order to account for the dynamics of the signal mode information, we first 
consider a general bipartite pure state of a system with fixed total energy that is 
composed of subsystems A and B, each with finite Hilbert space dimensions 6a and 
ds-, respectively. It is known that subsystem B will be nearly thermal and thus contain 
approximately no information as long as d& 3> ds, with the information content of 
subsystem B becoming apparent only after dA ~ g?b[IH IH2 HI]. Similar dynamics for 
the information content of the signal mode is shown in Fig. |4} However, the three 
harmonic oscillator modes from which a pure state of the trilinear Hamiltonian Eq. ((!]) 
is composed all contain an infinite set of states, preventing the direct application of 
dimensional arguments to our model. One can, however, define an effective subspace 
dimension for each model 



1 

d i ( r ) = Tr[of(T)] ' i = a,b,C (56) 
determined by the purity of the effective thermal state Cj(r) with temperature given by 
Eq. (53). This definition is motived by the fact that the purity of a mixed state density 
matrix is proportional to the number of basis states over which the fractional population 
of the mixed state is nonzero. For a state with an energy (quanta) constraint, using the 
thermal state cr(r) gives the minimal value for the purity, and as such Eq. (56) represents 
an effective maximum number of states constrained by the number of quanta in the mode 
at time r. Therefore, if we partition our initial pure state into bipartite subsystems d a 
and dbc composed of pump and combined signal and idler modes respectively, then the 
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initial effective subspace dimensions satisfy d<f ^> = (^ ff ) , where the equality is 
due to the symmetry between signal and idler modes. In Fig. |4]d we plot the times at 
which d e ^ = die f° r each of the initial pump mode coherent states considered in Fig. [4^,. 
Just as for finite dimensional pure states, the information contained in the signal or 
idler subsystems remains nearly zero until after rf„ ff (r) ~ G^f (r), provided we define the 



effective subspace dimension using Eq. (56). Our results are in agreement with those of 
Page|i 



where a similar argument was put forth for the evolution of information in the 
Hawking radiation from a finite dimensional black hole. 



4-4- Numerical Results for the Full Trilinear Evolution 

Unlike Hawking radiation, which is well-modeled using the quantum short-time 
approximation, backaction due to the build up of quanta in the signal and idler modes 



has a strong effect on the evolution of the trilinear Hamiltonian (16). In Fig. |2| we 
see that backreaction effects quickly lead to differing evolutions between the short-time 
and full quantum dynamics; backreaction prevents the full transfer of quanta from the 
pump mode[46j, and results in the oscillation of quanta between pump and signal/idler 
modes. To account for these backaction terms we repeat the analysis in Sec. |4.3| by 
numerically calculating the full dynamics of Eq. (17). Fig. [5] shows the modifications 
to both the fidelity and information content of the signal mode when backaction is 
included in the dynamics, as well as the subspace dimension condition d^(r) 



As expected, the short-time dynamics in Fig. [5] are nearly identical to those in Fig. [4] 
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Figure 5. a) Fidelities of the signal mode using numerical simulations of Eqs. (17) 
with a pump mode initially in a coherent state with occupation number (N a ) = 1 
(blue), 3 (red), 6 (green), and 9 (black), b) Corresponding information present in the 
signal mode spectrum and times at which <i° ff (r) = rf^(r) (dashed lines). 



up until d^(r) « dlc(r) which occurs later than in the short-time approximation as the 
backaction from particles in the signal mode begins to impede the transfer of quanta 
from the pump. The peaks in the information content of the signal mode correspond to 
the times at which backaction from the pump mode has completely stopped the transfer 
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of energy between modes. With the majority of quanta now in the signal and idler 
modes, the flow of energy reverses directions as the signal and idler now drive the pump 
mode leading to the oscillations seen in Fig. |5j As a model for black hole evaporation 
these oscillations represent the unphysical process of Hawking radiation flowing back 
into the black hole; the connection between the trilinear Hamiltonian and Hawking 
emission is valid only for the initial transfer of quanta from pump to signal/idler modes. 



5. Tripartite Entanglement 

With the pump mode quantized, one can consider the entanglement between pump and 
signal/idler modes. In Sec. [2j it was shown that ignoring this entanglement (pump 
mode treated classically) and tracing over one mode of a two-mode squeezed state leads 
directly to the thermal characteristics of the remaining system. In the full quantum 
description the statistics of the signal mode is obtained by tracing over both idler 
and pump modes and as such the distribution of entanglement between modes of 
this now tripartite system is important for characterizing the spectrum of the signal 
mode. Additionally, entanglement with the pump mode may appreciably alter the state 
dynamics of the pump compared to the classical approximations in Sees. [2] and [3j Given 
that multipartite entanglement is not as well understood as for the bipartite case, we 
begin by again considering the system to be partitioned into two subsystems consisting 
of the pump mode and the combined signal-idler. This bipartite separation into pump 
and signal-idler subsystems allows us to use the mutual information [42J 

la— be = S a + Sb c S a b c (57) 

as a measure of total correlations, both classical and quantum, between subsystems . 
Since our tripartite system state remains pure throughout its evolution, the total entropy 
S a bc = and the subsystem entropies satisfy S a = Sb c - The mutual information is 
therefore twice the entropy of the pump mode subsystem, I a -bc = 2S a . Thus, the 
subsystem entropy of the pump is a direct measure of entanglement with the signal and 
idler modes, which is not taken into account in either the parametric or semi-classical 
solutions. Of course the signal and idler subsystems are also entangled with each other 
as in Sees. [2]and[3j With both signal and pump modes in identical states, the signal-idler 
mutual information is given by 

h-c = S + S c — Sbc — 2Sb — S a (58) 

indicating the tradeoff between the entanglement of the signal/idler and that of the 
pump mode Eq. (57). In Fig. [6^ we plot the mutual information Eqs. (57) and (58) for 



a pump initially in a coherent state of amplitude (N a (0)} = 9 for the parametric, semi- 
classical and short-time approximations along with the full numerical solution obtained 



using Eqs. (17). We see the mutual information (pump entanglement) l a -bc begins to 
become appreciable around the same time as the information content of the signal mode 
becomes apparent in Figs. |4Jd and|5)D indicating the increasing role of the pump mode in 
the dynamics. The numerical solution shows that this increase in pump entanglement 
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reduces the signal-idler mutual information Jfe_ c with respect to the semi-classical and 



parametric approximations as expected from Eq. (58). In the full quantum dynamics 



the pump is never allowed to be fully depleted and is therefore always entangled with 
the signal and idler. However in the short-time approximation the late-time pump mode 
is depleted and approaches the ground state where S a = 0. 

Finally, we show that the entanglement buildup with the pump mode given by 
I a -bc not only affects the signal and idler modes, but also results in amplitude dependent 
squeezing of the pump. In Fig. we plot the amplitude of the squeezing parameters [M]. 

a(r) ± a + (r) 



Q± 



AXl 



1, 



(59) 



for the pump mode, where ^AX±j is the variance of the quadratures. As can be seen, the 
entanglement with the signal and idler modes leaves the pump mode in a non-classical 
squeezed state over the time range of interest. 
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Figure 6. a) Signal-idler mutual information I\,-c f° r pump with (iV a (0)) = 9 
in the parametric (red), semi-classical (blue), quantum short-time (green), and full 
numerical solution (black). Dashed lines represent mutual information I a -bc indicating 
entanglement with the pump mode, b) Quadrature squeezing parameters q± of the 
pump mode given the pump is initially in a coherent state with occupation number 
(N a ) = 1 (blue), 3 (red), 6 (green), and 9 (black). Negative values indicating squeezing 
of the q- quadrature. Dashed lines represent the corresponding q + quadrature. 



6. Conclusion 

In this paper we have investigated the effect of a dynamical quantized pump mode on the 
generation of quanta in a parametric amplifier. We have shown that a quantized pump 
mode leads naturally to a non-thermal spectrum for the signal and idler modes once 
the pump has released nearly half of its initial energy, such that the effective subspace 
dimensions of the pump and signal/idler mode systems approximately coincide. The 
departure from a thermal state indicates that the signal spectrum contains information 
that may be used to partially reconstruct the initial state of the pump mode. Once 
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quantized, the pump mode becomes entangled with the signal and idler leading to 
non-classical squeezed states of the pump. As a simple, zero-dimensional model of 
the Hawking effect, the present findings lend support to the possibility for non-thermal 
emission from a quantum (as opposed to semiclassical) black hole; the emitted radiation 
is entangled with the quantized gravitational degrees of freedom and yields information 
concerning initial formation of the black hole. 
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